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(54) A method for synchronizing a receiver, a system, and an electronic device 



(57) The invention relates to a method for synchro- 
nizing a receiver (1) with a transmitted code-modulated 
spread spectrum signal, tn the method, at least one ref- 
erence code (r(x)) is used, which corresponds to a code 
used in the modulation, and the frequency shift of the 
transmitted signal and the code phase of the code used 
in the modulation are determined. In the method, a cor- 
relation step is taken to form a correlation function matrix 
on the basis of the received signal and said reference 



code. The step is reiterated a plurality of times. The it- 
erations results are added non-coherently to form a non- 
coherent search matrix. Elements in said non-coherent 
search matrix are modified on the basis of at least one 
statistical property of the elements in said non-coherent 
search matrix, and/or high-pass filtering is performed 
before the formation of said correlation function matrix. 
The invention also relates to a system and an electronic 
device in which the method is applied. 
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Description 

[0001] The present invention relates to a method for synchronizing a receiver with a transmitted code-modulated 
spread spectrum signal, in which method at least one reference code is used, which corresponds to a code used in 

5 the modulation, and the frequency shift of the transmitted signal and the code phase of the code used in the modulation 
are determined. The invention also relates to a system comprising a receiver with acquisition means for acquisition of 
a transmitted code-modulated spread spectrum signal, means for generating at least one reference code corresponding 
to the code used in the modulation, and means for determining the frequency shift of the transmitted signal and the 
code phase of the code used in the modulation. The invention further relates to an electronic device comprising a 

10 receiver with acquisition means for acquisition of a transmitted code-modulated spread spectrum signal, means for 
generating at least one reference code corresponding to the code used in the modulation, and means for determining 
the frequency shift of the transmitted signal and the code phase of the code used in the modulation. 
[0002] Particularly in the reception of satellite signals, a significant problem is that the signal strength is low in the 
receiver. Thus, noise and other interference may significantly disturb the signal reception and, in some situations, it 

15 may be even impossible to receive signals. In addition to noise, interference may be caused by other signals whose 
carrier frequency is the same or almost the same as the frequency of the signal to be received. Furthermore, a DC 
voltage component in the signal may cause interference in direct conversion receivers and receivers applying a low 
intermediate f reque* y A spurious signal at the receiving frequency and the DC shift may be shifted and/or aliased at 
different steps of th> processing to spurious signals at different frequencies. It is difficult to eliminate these different 

20 sources of interference mentioned above, and moreover, the elimination of different forms of interference is not usually 
successful by applying a single interference suppression method. 

[0003] There are known systems based on the code division multiple access (CDMA). In such systems, the infor- 
mation to be transmitted is modulated with an individual code for each transmitter, a so-called pseudo random se- 
quence. Thus, the same frequency can be used as the carrier frequency in different transmitters. As a result of the 

25 modulation, a code-modulated wideband signal is generated. This signal is received in the receiver, and an attempt is 
made to synchronize the receiver with it. The receiver knows the code used in the transmitter and uses it in the acqui- 
sition of the signal. This code can thus be used to distinguish between signals from different transmitters, even though 
the carrier frequencies were substantially the same. In the acquisition, a correlation technique is normally used to 
correlate the received signal with a code corresponding to the code used by the transmitter and generated in the 

30 receiver. The correlation result is examined to find the timing and/or frequency of the incoming signal corresponding 
to maximum correlation, i.e. the best alignment. However, interference may distort the correlation result or cause false 
maximum points, wherein the acquisition is not necessarily successful. For this, solutions have been developed to find 
the correct timing and frequency e.g. by prolonging the correlation time. However, particularly with weak signals, such 
as when signals transmitted by satellites are received indoors, acquisition by a conventional receiver would require 

35 several hours or even days, which, in practical situations, would make it even impossible to use such receivers indoors. 
[0004] One known system applying the CDMA technology is the Global Positioning System (GPS) comprising several 
satellites orbiting the earth. Each operating satellite of the GPS system transmits a so-called L1 signal at the carrier 
frequency of 1575.42 MHz. This frequency is also indicated with 154f 0 , where f 0 = 10.23 MHz. Furthermore, the satellites 
transmit another ranging signal at a carrier frequency of 1 227.6 MHz called L2, i.e. 120f 0 . In the satellite, these signals 

40 are modulated with at least one pseudo sequence. This pseudo random sequence is different for each satellite. In 
each satellite, for modulating the L1 signal, the pseudo random sequence used is e.g. a so-called C/A code (Coarse/ 
Acquisition code), which is a code from the family of the Gold codes. Each GPS satellite transmits a signal by using 
an individual C/A code. The codes are formed as a modulo-2 sum of two 1023-bit binary sequences. The first binary 
sequence G1 is formed with the polynomial X 10 + X 3 + 1 , and the second binary sequence G2 is formed by delaying 

45 the polynomial X 10 + X 9 + X 8 + X 6 + X 3 + X 2 + 1 in such a way that the delay is different for each satellite. This 
arrangement makes it possible to generate different C/A codes by using identical code generators. The C/A codes are 
thus binary codes whose chipping rate in the GPS system is 1.023 Mchips/s. The C/A code comprises 1023 chips, 
wherein the repetition interval (epoch) of the code is 1 ms. The carrier of the L1 signal is further modulated by navigation 
information at a bit rate of 50 bit/s. The navigation information comprises information about the "health", orbit, time 

50 data of the satellite, etc. In the GPS system, the codes used in the modulation of the L1 signal are not particularly 
efficient in view of eliminating the above-presented narrow-band interference. Thus, the cross-correlation caused by 
a strong spurious signal in the weaker signal to be received may prevent the receiver from acquiring this signal to be 
received. 

[0005] The receiver must perform the acquisition e.g. when the receiver is turned on and also in a situation in which 
55 the receiver has not been capable of receiving the signal of any satellite for a long time. Such a situation may easily 
occur e.g. in portable devices, because the device is moving and the antenna of the device is not always in an optimal 
position in relation to the satellites, which impairs the strength of the signal coming in the receiver. Also, in urban areas, 
buildings affect the signal to be received, and furthermore, so-called multipath propagation may occur, wherein the 
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transmitted signal comes into the receiver along different paths, e.g. directly from the satellite (line-of-sight) and also 
reflected from buildings. Due to this multipath propagation, the same signal is received as several signals with different 

phases. 

[0006] The above-mentioned acquisition and frequency control process must be iterated for each signal of a satellite 
5 received in the receiver. Consequently, this process takes a lot of time, particularly in a situation, in which the signals 
to be received are weak. To speed up this process, some prior art receivers use several correlators, wherein it is 
possible to search for several correlation peaks simultaneously. In practical solutions, the process of acquisition and 
frequency control cannot be accelerated very much solely by increasing the number of correlators, because the number 
of correlators cannot be increased infinitely. 
10 [0007] In some prior art GPS receivers, FFT technique has been used in connection with conventional correlators 
to determine the Doppler shift of the received GPS signal. These receivers use the correlation to restrict the bandwidth 
of the received signal to 1 kHz. This narrow-band signal is analyzed with FFT algorithms to determine the carrier 
frequency. 

[0008] It is an aim of the present invention to provide a method for accelerating the acquisition, and a receiver in 
*5 which the acquisition of the transmitted signal can be performed significantly faster than in receivers of prior art, also 
with weak signals and interference. The invention is particularly suitable for use in positioning receivers but also in 
other receivers, preferably CDMA receivers, in which the receiver must be synchronized with a spread spectrum signal. 
The invention is based on the idea that a non-coherent search matrix formed of a correlation function matrix is post- 
processed by statistical methods. Furthermore, it is possible to carry out high-pass filtering of the computed values of 
20 the correlation functions for the output signal. To put it more precisely, the method according to the present invention 
is primarily characterized in that a correlation step is taken to form a correlation function matrix on the basis of the 
received signal and said reference code, a non-coherent search matrix is formed of said correlation function matrix, 
and that the method comprises at least one of the following steps: 

25 - modifying the elements of said non-coherent search matrix on the basis of at least one statistical property of the 
elements of said non-coherent search matrix, 

performing high-pass filtering before forming said correlation function matrix. 

[0009] The system according to the present invention is primarily characterized in that the system also comprises 
30 correlation means for forming a correlation function matrix on the basis of the received signal and said reference code, 
means for forming a non-coherent search matrix from said correlation function matrix, and at least one of the following 

means: 

means for modifying the elements of said non-coherent search matrix on the basis of at least one statistical property 
35 of the elements of said non-coherent search matrix, 

means for performing high-pass filtering before the formation of said correlation function matrix. 

[001 0] Further, the electronic device according to the present invention is primarily characterized in that the electronic 
device also comprises correlation means for forming a correlation function matrix on the basis of the received signal 
*o and said reference code, means for forming a non-coherent search matrix from said correlation function matrix, and 
at least one of the following means: 

means for modifying the elements of said non-coherent search matrix on the basis of at least one statistical property 
of the elements of said non-coherent search matrix, 
45 - means for performing high-pass filtering before the formation of said correlation function matrix. 

[001 1 ] Considerable advantages are achieved by the present invention when compared with methods, systems and 
electronic devices of prior art. By the method according to the invention, it is possible to significantly suppress the 
interference to DC voltage shift (offset) or due to other strong in-band signals. It is thus possible to receive also weaker 
50 signals than by receivers of prior art. 

[0012] In the following, the present invention will be described in more detail with reference to the appended drawings, 
in which 



55 



Fig. 1 shows, in a simplified block chart, a receiver in which the method according to the invention can be applied, 

Fig. 2 shows the correlation step of the method according to a preferred embodiment of the invention in a reduced 
chart, 
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10 



Fig. 3 shows the analysis step of the method according to a preferred embodiment of the invention in a reduced chart, 

Fig. 4 shows the summing step of the method according to a preferred embodiment of the invention in a reduced 
manner, 

Fig. 5 shows the conclusion step of the method according to a preferred embodiment of the invention in a reduced 
manner, and 

Fig. 6 shows an electronic device according to a preferred embodiment of the invention in a reduced block chart. 



[0013] In the following description, the method according to an advantageous embodiment of the invention will be 
exemplified with a receiver shown in Fig. 1 and intended for use as a receiver of signals transmitted by GPS satellites, 
for positioning. However, it will be obvious that the invention is not limited solely to such receivers, but the receiver 
used can also be another receiver intended for receiving a code modulated signal. Also, the receiver may differ from 

15 the receiver of Fig. 1 in detail. 

[0014] In the receiver 1 of Fig. 1 , the signal to be received is preferably converted into an intermediate frequency in 
a converter block 2. At this stage, the signal comprises two components, known as such: the I and Q components, 
with a phase difference of 90° therebetween. These analog signal components converted to an intermediate frequency 
are digitized in a digitization block 3 and led to a multiplier block 4. In the multiplier block 4, the I and Q components 

20 of the digitized signal are multiplied by a signal generated by a numerically controlled oscillator (NCO) 5. This signal 
of the numerically controlled oscillator 5 is intended to correct the frequency offset caused by the Doppler shift and the 
frequency error of the local oscillator (not shown) of the receiver t . The signal generated by the multiplier block 4 is 
led to an acquisition block 6. This acquisition block 6 finds the spreading code phase and frequency offset of the 
satellite, to be used in operations after the acquisition. During the acquisition, a control block 7 is used to control a 

25 scanning block 8, by means of which the frequency of the numerically controlled oscillator 5 is adjusted, if necessary. 
The control block 7 controls a first switch 9 to switch either the signal generated by this scanning block 8 to the nu- 
merically controlled oscillator 5 during the acquisition, or the control signal generated by a code determination block 
11 to the numerically controlled oscillator 5 after the acquisition. A second switch 10 is used to control this operation 
of the tracking block 11. This tracking block 11 forms a part of the code phase locked loop, known as such, and the 

30 carrier phase locked loop (not shown). 

[0015] After turning on the operating voltages, or in a situation in which the receiver 1 has not been capable of 
receiving the signal of GPS satellites for a long time, the receiver 1 performs a two-dimensional searching step for 
each satellite whose signal is to be received. In this two-dimensional searching step, the aim is to determine the carrier 
frequency and code phase for each satellite. The carrier frequency is thus affected by the Doppler shift caused by the 

35 movement of the satellite, as well as by inaccuracies in the local oscillator of the receiver. The frequency uncertainty 
may be relatively large, for example ± 6 kHz, wherein the receiver 1 must scan a frequency range of about 12 kHz in 
relation to the actual transmission frequency (L1 = 1575.42 MHz). Moreover, the receiver 1 does not know the precise 
code phase, wherein the receiver must also determine the code phase from 1023 possible code phases. This will result 
in a two-dimensional search process, in which one dimension is the frequency offset in the range of 12kHz and the 

40 second dimension is the code phase from 1 023 different code phases. In the method according to a preferred embod- 
iment of the invention, it is possible to scan a frequency range of about 500 Hz at a time., wherein the method is iterated, 
if necessary, 24 times to cover the whole frequency range of 12 kHz to be scanned. It is obvious that the example 
values in this description are only used to clarify the invention, but not as restricting examples. The invention can also 
be applied in systems other than GPS systems, wherein e.g. said frequency values, code phases and the number of 

45 codes may vary. The search process does not necessarily need to be two-dimensional, but in some embodiments, it 
may also have one dimension or more than two dimensions. Thus, in a single-phase search process, the aim is to find 
out, for example, the code phase. 

[0016] We shall now describe the acquisition operation based on forming a two-dimensional search matrix in the 
receiver 1 of Fig. 1. To start the acquisition, the scanning block 8 sets the frequency of the numerically controlled 

50 oscillator 5 so that the receiver 1 preferably receives the lowest frequencies in the frequency range, in this example 
from 1575.414 MHz to 1575.4145 MHz. The receiver may also determine the starting frequency so that the receiver 
utilizes e.g. previously determined position data and/or almanac data, wherein the acquisition can be further acceler- 
ated. Samples of the received signal are preferably stored in a sample vector forming block 12 as complex sample 
vectors p k (1). p K (2)...p k (N) ; which each comprise 1023 samples in this preferred embodiment. In this advantageous 

55 embodiment, the frequency of storing samples in the sample vector forming block 12 is substantially the same as or 
higher than the chipping rate : that is, 1.023 to 10.23 Msamples/s. The sample vectors are immediately adjacent so 
that the next sample vector immediately follows in time the previous sample vector, i.e. the time difference between 
the last sample of a sample vector and the first sample of the next sample vector is substantially the same as the time 
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difference between successive samples in the sample vector. Thus, these 1023 samples represent a signal with the 
length of about 1 ms, wherein in the time-to-frequency conversion, the frequency range is about 1 kHz, of which a part 
can be utilized. The sample vector forming step is indicated by reference 101 in the appended Fig. 2. 
[0017] The number of sample vectors is advantageously N, in which N is preferably a power of two. Furthermore, in 

5 an advantageous embodiment of the invention, the forming of sample vectors is iterated K times, as will be presented 
below in this description. Subindex k will be used hereinbelow to indicate different iterations. When determining the 
value of the number N of the sample vectors p k (1), p k (2)...p k (N) in the GPS system, one must take into account that 
the signal contains information modulated at a bit rate of 50 bits/s by binary phase modulation. Another factor limiting 
this number N of the sample vectors p k (1 ), p k (2)...p k (N) is the frequency stability of the local oscillator in the receiver. 

w [0018] In addition to the step of forming the sample vectors, in the acquisition method according to the invention, a 
correlation step is taken in order to form a correlation function matrix. This correlation step can be taken partly already 
during the sampling, or after the formation of N number of sample vectors p k (1), p k (2)...p k (N). If the correlation step is 
taken e.g. so that, after the storage of each sample vector, a time-to-frequency conversion is computed for it, such as 
a Fast Fourier transform (FFT), the same time-to-frequency converter can be used for all N sample vectors p k (1), p k 

15 (2)...p k (N). However if the correlation step is taken after the storage of the N sample vectors, either a separate time- 
to-frequency converter must be used for each sample vector, or the time-to-frequency conversions are performed for 
the different sample vectors one after the other in the same time-to-frequency converter. Figure 2 shows the correlation 
step of the method, wherein a correlation function matrix C x k is formed of the sample vectors p k (1 ), p k (2).. .p k (N). Below 
in this description, the Fourier transform and the inverse Fourier transform will be primarily used as examples of the 

20 time-to-frequency conversion and the inverse conversion, i.e. frequency-to-time conversion, respectively; however, it 
will be obvious that the present invention is not limited solely to these examples. 

[0019] Each sample vector p k (1), p k (2)...p k (N) is subjected to discrete Fourier transform 102 ; preferably fast Fourier 
transform (FFT): 
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P k (i)=FFT(p k (i)), (1) 



in which i = 1 N 

[0020] This is illustrated by blocks FFT1 , FFT2, .„, FFTN in Fig. 2. In arithmetic operations in practice, preferably 
1 024 values are used, because the discrete Fourier transform can thus be implemented in a considerably more efficient 
way in practical applications (with the FFT algorithm) than when using 1023 values. One way of achieving this is to 
add an extra zero to make the number of elements 1 024. This will have an insignificant effect on the transform result. 
[0021] The codes of all those transmissions, for whose reception the receiver 1 is intended, are preferably stored in 
the receiver. Thus, in connection with the GPS system, the receiver 1 contains the stored reference codes r(x) corre- 
sponding to the C/A codes for the satellites (not shown) in the GPS system, wherein x refers to the satellite index, x 
ranging for example from 1 to 36. The reference codes do not necessarily need to be stored but they can also be 
generated in the receiver. At the correlation stage, the reference code of the satellite which transmitted the signal to 
be acquired by the receiver is selected or generated. The reference code is reversed in time. This reversed reference 
code : which is indicated with reference r(x) in Fig. 2, is subjected to a discrete Fourier transform 103, preferably a Fast 
40 Fourier transform, or FFT: 

R(x)=FFT(7(x)) (2) 



[0022] The time-reversed reference code r(x) and/or its fast Fourier transform may already have been stored in 
advance in the memory means of the receiver, or it is formed from the reference code r(x) in connection with the 
acquisition. 

[0023] Next, in the correlation step, the Fourier transform result P k (i) of each sample vector p k (i) is multiplied 104 
with the Fourier transform r(x) result of the inverse reference code R(x) : 

Mx.k(0 = R(x)-P k (i) (3) 

[0024] These products of the multiplications are further subjected to an inverse Fourier transform 1 05, wherein the 
result is the cross correlation of the reference code r(x) and the received signal with all possible integer delays (1023). 
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m xk (i)=iFFT(M xk (i)) (4) 

[0025] This is based on the fact that the Fourier transform of the convolution of the signals in the time domain cor- 
5 responds to the multiplication of the Fourier transformed signals, i.e. signals of the time domain converted to the fre- 
quency domain. As also the time-reversed reference code is used, the Fourier transform can be used to perform a fast 
correlation in discrete time. Thus, in this preferred example, the cross correlation result comprises 1023 elements. 
These cross-correlation results m x k (i) formed of different sample vectors p k (i) are compiled to a correlation function 
matrix C x k , in which the number of rows is the number N of the sample vectors. 
10 [0026] It is obvious that, instead of time-reversal of the reference code in time, time-reversed sample vectors p k (i) 
can be formed of the sample vectors p k (i), wherein the reference code r(x) is used directly and time-reversed sample 
vectors are used in the above-presented arithmetic operations. 

[0027] It should also be pointed out that in view of applying the present invention as a whole, the method used for 
producing the cross-correlation result is not significant as such. 
15 [0028] The rows in the correlation function matrix C x k formed in the correlation step represent the cross-correlation 
of the received signal and the reference code with various phase differences taken at intervals of one millisecond. The 
correlation function matrix can thus be presented as the following formula: 



20 



25 



35 



m x *(2) 
m Xtk (/V) 



(5) 



[0029] In the next step, i.e. the analysis step, a transposition 1 06 of the correlation function matrix C x k will be used, 
the rows representing samples of the signal in the time domain. Each row corresponds to a specific code phase dif- 
30 ference between the received signal and the reference code. Each row of the transpose of this correlation function 
matrix C x k is subjected to a Fourier transform 107 to form a coherent search matrix A x k . This is illustrated in the 
appended Fig. 3. 



A x . k =FFT( C y (6) 



40 [0030] In practical applications, a separate transposed matrix does not need to be formed from the correlation function 
matrix, but the elements of the stored correlation function matrix C x k are read from a memory 16 (Fig. 6) in a different 
direction, preferably in columns. 

[0031] The correlation function matrix C x k can also be formed e.g. by using matched filters, known as such. Thus, 
the output signal of the matched filter is preferably sampled by forming a sample signal suitable for the further processing 
45 steps, such as the Fourier transform 1 07. In this step, the sampling frequency is about 1 kHz, in the case of a receiver 
intended to receive signals from GPS satellites. 

[0032] Furthermore, the acquisition block 6 comprises a non-coherent summing step to improve the signal-to-noise 
ratio. To implement the non-coherent summing step, the above-presented steps of sample vector formation, correlation 
and analysis are iterated 108 K times (Fig. 4). This number K of iterations is preferably selected so that the signal-to- 

50 noise ratio can be improved to a sufficient degree but within a reasonable time. At each time of performing the analysis, 
one coherent search matrix \ k is formed, for which the non-coherent summing is performed to form a non-coherent 
search matrix S x . The non-coherent search matrix S x is preferably formed in the following way. From the complex 
elements a x k {ij) of each coherent search matrix A x k , preferably the magnitude (absolute value) or another numerical 
value is calculated, such as the second power of the magnitude of the element. The numerical values calculated from 

55 the corresponding elements of the non-coherent search matrix are summed up 1 09, i.e. the matrices are summed up, 
which can be presented by the formula: 
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]a M (1023,1]| ... |a„,,(l023,A/)| 



(7) 



10 [0033] In practical solutions, the non-coherent search matrix can be formed in at least two ways. In the first imple- 
mentation alternative, the coherent search matrix formed at each iteration time is stored. After the necessary iterations, 
the non-coherent search matrix is formed by summing up the corresponding elements according to Formula 8. In this 
implementation alternative, a memory is needed for storing ail the elements of the coherent search matrices. According 
to the second implementation alternative, one coherent search matrix is calculated first, and its values are copied as 

*5 elements of the sum matrix. Each following iteration time, a coherent search matrix is formed, whose values are added 
to the corresponding elements in the non-coherent search matrix. In this alternative, the summation of the correspond- 
ing elements is thus performed each iteration time. Thus, only one coherent search matrix is stored, wherein less 
memory space is required than in the first alternative. 

[0034] After the necessary iterations have been completed, the values of the elements s x (i,j) of this non-coherent 
20 search matrix S x are modified in a modification step in order to improve the rate and accuracy of the acquisition. In 
connection with an advantageous embodiment of the invention, this is implemented in the following way: The non- 
coherent search matrix is preferably used to form two filtration matrices to be used as the basis for the acquisition. To 
compute the filtration matrices, preferably the first and the second scaling factor are computed from the elements of 
the different code phases corresponding to each frequency, by using suitable statistical methods. This computation is 
thus performed according to columns in the case that each row in the search matrix corresponds to a given code phase 
difference between the received signal and the reference code. The first scaling factor is computed, for exampie, as 
the average of all the different code phases corresponding to one frequency. Next, this first scaling factor is subtracted 
from all the values in this column, and these results are compiled to one column in the first filtration matrix. The second 
scaling factor can be preferably computed by means of the standard deviation. Thus : the standard deviation is calcu- 
30 lated from the values in each column and can be scaled by multiplying it with a suitable number, for example seven. 
Next, this second scaling factor is subtracted from the column values and the resulting values are stored as column 
values of the second filtration matrix. If any value is negative after said subtractions, it is preferably set to zero. 
[0035] The above-presented steps are iterated for each column of the non-coherent search matrix; that is, separate 
first and second scaling factors are computed for each column to be used for modifying the values of the column in 
35 question, and are stored as the column values for the first and second filtration matrices. After forming the first and 
second filtration matrices, the aim is to find from the first and/or second filtration matrices a value which exceeds a 
predetermined threshold value and is clearly greater than the other values. If such a value is found, it indicates the 
code phase difference as well as the frequency offset, because it is probably a signal transmitted by a satellite. If the 
signal is not the desired signal but noise or another spurious signal, no significant correlation peaks should occur. The 
40 code phase difference and the frequency difference are manifested by the row index and the column index of this 
highest value, respectively. The above-presented modification based on the average will efficiently eliminate DC shift 
interference. This is based on the fact that the DC shift will have a substantially equal effect on all the values in a 
column. Thus, the average of the values in such a column is relatively high. In a corresponding manner, the modification 
based on the standard deviation will effectively affect cross-correlation interference, such as cross-correlation results 
45 caused by a strong spread spectrum transmission at substantially the same frequency. This is based on the fact that 
cross-correlation caused by a spurious signal will only affect some of the different code phases, wherein the standard 
deviation and also the variance will be relatively high. The above-mentioned computations can be preferably made in 
a digital signal processor 15 or in another computing means suitable for the purpose. 

[0036] The example of Fig. 5 shows a situation in which interference is suppressed and a value clearly higher than 
the other values has been found. However, if no such value is found in the non-coherent search matrix s x, i.e. a signal 
transmitted by the searched satellite was probably not received in the frequency range under examination, the fre- 
quency band under examination is changed and the above-presented steps are taken to form a non-coherent search 
matrix. By this method, the whole range of 1 2 kHz under examination can be scanned by iterating the above-presented 
steps a sufficient number of times. 

[0037] If necessary, the above-presented steps can be iterated for the whole frequency band under examination, 
storing the filtration matrices formed at the different iteration times, or only the possible peaks, before searching for 
the highest correlation peak. In this way : the possibility of misinterpretations can be reduced e.g. in such a situation in 
which the threshold value is set too low and a spurious signal not yet filtered enough may cause a misinterpretation. 
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[0038] Naturally, the invention can also be applied by forming one filtration matrix first, for example the first filtration 
matrix. If one value clearly higher than the other values is not found from this filtration matrix, the second filtration 
matrix will also be formed and the search will be carried out on the basis of its values. It is also possible to form one 
filtration matrix first, for example the first filtration matrix, and if no value substantially higher than the other values is 
5 found from this filtration matrix, the second filtration matrix can be formed of this first filtration matrix, to carry out the 
search on the basis of its values. 

[0039] In the method according to another advantageous embodiment of this invention, the high-pass filtering 110 
(Fig. 2) is performed before the formation of the search matrices. This high-pass filtering is preferably performed after 
the computation of the correlation function, before the formation of the correlation function matrix C x k . Thus, the high- 
10 pass filtering is performed for the cross -correlation results m x k (i) formed of the sample vectors p k (i), and the correlation 
function matrix C x k is formed of the high-pass filtered values. 

[0040] We shall now discuss the effect of high-pass filtering on the signal. The DC shift in the output of the analog- 
to-digital converter 3 is converted to a low-frequency sinusoidal signal in a mixer 4. The frequency of this sine wave 
to be formed is substantially the same as the frequency of the numerically controlled oscillator 5. Assuming that the 

15 frequency shift at the carrier frequency is limited to be less than 100 kHz, the frequency of the sine wave will always 
be lower than 1 00 kHz. The desired signal in the output of the means for computing the correlation function has sub- 
stantially the length of one scatter code symbol (chip) which, in the GPS system, thus means 0.98 us impulses (= 
1/1 023000 s) This impulse contains considerably higher frequencies than the sine waves generated from the DC shift. 
Thus, the high-pass filtering will not suppress the desired signal to an adverse extent. 

20 [0041] In this context, it should be mentioned that low-frequency interference will always act as said DC shift. Thus, 
the effect of any interference at a frequency lower than the barrier frequency of the high-pass filter can be eliminated 
by the high-pass filter. Said position of the high-pass filter after the computation of the correlation function is advanta- 
geous, because the word length of the samples there is normally at least ten bits, which will be sufficient to prevent 
the high-pass filter from impairing the quality of the actual signal. However, the resolution of the analog-to-digrtal con- 

25 verter 3 is, e.g. for reasons of power consumption and to limit the complication of the required equipment, normally 
only two to four bits, which is normally not an adequate word length to perform high-pass filtering of sufficient quality. 
[0042] In addition to the high-pass filtering, if necessary, one filtration matrix may be advantageously formed of the 
search matrix, preferably using the standard deviation as the statistical function, as presented above. In this way, it is 
possible to suppress both DC shift interference and cross-correlation interference substantially more efficiently than 

30 in receivers of prior art. 

[0043] After determining the correct frequency offset and code phase, the receiver can be set in a tracking mode. 
With the weakest signals, data reception will not be successful, but, in a way known as such, one must turn to data 
obtained e.g. via a mobile communication network. Distance measuring is still possible with the reduced accuracy. The 
tracking mode is set by changing the switches (Fig. 1) to another position, wherein the received information is led to 

35 the tracking block 11 , in which also a feedback is formed for the fine adjustment of the frequency of the numerically 
controlled oscitlator5. 

[0044] To calculate the position, the receiver performs the reception of a signal preferably on the basis of the signal 
received from at least four satellites. Thus, the above-presented acquisition is iterated, if necessary, for the signal of 
each satellite, wherein the reference sequence r(x) is selected to be the code of the satellite to which the acquisition 
40 is made at the time. 

[0045] Furthermore, we shall discuss the effect of aliasing in the receiver 1 according to the invention. The aliasing 
phenomenon occurs, if signals whose frequency is higher than one half of the sampling rate are led to the analog-to- 
digital conversion. As mentioned above, the DC shift in the output of the analog-to-digital converter 3 is converted to 
a low-frequency sinusoidal signal by the frequency conversion in the mixer 4. The frequency of this sinusoidal signal 

45 is substantially the same as the negation of the frequency of the numerically controlled oscillator 5. For example, if the 
frequency shift is +20 kHz and the search is being carried out at a frequency range with a medium frequency of -1300 
Hz, the frequency of the numerically controlled oscillator 5 is set to the value +20000 + -1300 Hz = +18700 Hz. The 
frequency conversion converts the frequency of the satellite signal to be received to about 0 Hz, wherein the frequency 
of the sine wave caused by the DC shift is shifted downwards to about -1 8700 Hz. After the frequency conversion and 

50 the decimation, the signal is subjected to the computation of the correlation function. The computation of the correlation 
function suppresses signals of different frequencies in different ways, wherein the effect of the DC shift interference 
depends on the magnitude of the frequency shift. 

[0046] After the computation of the correlation function, the signal is sampled again at a sampling frequency of, for 
example, 1 kHz. At this stage, all the signals outside the frqeuency of +/- 500 Hz will be aliased to the range from -500 
55 Hz to +500 Hz. Also, the signal formed of the DC shift will be aliased to this range. The frequency of this aliased signal 
formed by the DC shift can be computed by subtracting or adding a suitable quantity of whole kilohertzes to the fre- 
quency of the sinusoidal signal. For example, the frequency -18700 of the sinusoidal signal will be 300 Hz after the 
aliasing, because -18700 + 19000 = +300. This can be presented by the formula W^ong' 1000 *! fork/ 1000 !. in which 
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20 



the denotation Q refers to rounding to the closest Integer, f orig = -f NCO , and f NCO is the frequency of the numerically 
controlled oscillator 5 in the frequency conversion. Using this formula, it is possible to compute the frequency of the 
interference caused by the DC shift, wherein it is possible to avoid the column corresponding to this frequency in the 
non-coherent search matrix. 

5 [0047] Also, the effect of cross-correlation can be analyzed in a corresponding manner by means of the aliasing 
phenomenon. The frequency conversion will also cause shifting and aliasing of the frequency of the cross-correlation 
interference. For example, if the value of the frequency shift is +20 kHz, the Doppler frequency of the satellite causing 
cross-correlation is +2400 Hz, and the medium frequency of the frequency range under examination is -1300 Hz, the 
frequency of the numerically controlled oscillator 5 is set to the value 18700 Hz. Thus, the frequency of the cross- 

10 correlation signal will be changed to +3700 Hz. 

[0048] Due to the properties of the Gold codes, the computation of the correlation function will amplify not only the 
desired signal but also the cross-correlation signal. However, the amplification of cross-correlation signals is typically 
about 20 to 25 dB lower than the amplification of the desired signal. Thus, the cross-correlation will primarily cause 
problems when the signal strength of the interfering satellite is more than 20 dB higher than the strength of the desired 

*5 signal. 

[0049] The aliasing will also affect the cross-correlation . By resampling at the frequency of 1 kHz, the cross-correlation 
signals will be aliased to the range from -500 Hz to +500 Hz. The above-presented formula for calculating the aliasing 
frequency will also apply in this case. In the formula, the original signal f orig is now represented by the frequency of the 
cross-correlation signal after the frequency conversion. 

[0050] A majority of the blocks required for implementing the method can be implemented, for example, in a digital 
signal processor (DSP) . When making the Fast Fourier Transforms, it is possible to use either hardware-based solutions 
or software implementations of the digital signal processor. Furthermore, to control the operation of the receiver, it is 
possible to use a control means, preferably a microprocessor or the like. In the high-pass filtering, a hardware-based 
implementation is preferably used, because the computation time available is short. For example, in a GPS application, 
25 the sampling frequency in the output of the means for computing the correlation function is in the order of 1 MHz to 10 
MHz. 

[0051] Further, the appended Fig. 6 shows an electronic device 24 according to an advantageous embodiment of 
the invention, comprising a wireless communication device and a positioning receiver. A first antenna 13 is used to 
receive a signal transmitted from positioning satellites. The received signal is transferred to a first radio part 14, in 
which the signal is converted to an intermediate frequency and digitized. A first radio part comprises, for example, a 
converter block 2, a digitizing block 3 and a multiplication block 4 for the receiver of Fig. 1 . The digitized signal, which 
in this case preferably comprises I and Q components, is led to a digital signal processor 15 for e.g. forming sample 
vectors in a sample vector formation block 12. Samples are stored in first memory means 16 which comprise e.g. a 
random access memory and preferably also a read-only memory and/or a non-volatile random access memory for the 
storage of the program code of the digital signal processor 1 5. In this embodiment, e.g. operations of the acquisition 
block 6, such as the formation of the correlation function matrix C x k e.g. with Fourier converters FFT1, FFT2, .... FFTN 
and/or matched filters, are implemented in the signal processor 15. Preferably, the coherent search matrix A x k is formed 
and the non-coherent summing step and the modification step are also taken in the digital signal processor 15. The 
digital signal processor 15 transmits information about the computed phase difference and the frequency offset to the 
processor block 1 7 which comprises e.g. a microprocessor and I/O logic. The processor block 1 7 controls the scanning 
block 8 as well as the first switch 9. Preferably, the tracking block 1 1 is at least partly impleme nted as program commands 
of the processor block. Second memory means 18 are used as the data memory and the program memory for the 
processor block 17. It is obvious that the first memory means 1 6 and the second memory means 18 can also comprise 
a common memory. Positioning information can be displayed to the user on a display 19. In the electronic device of 
Fig. 6, the high-pass filtering is implemented on hardware basis in a high-pass filter 110. 

[0052] Also, functions of a wireless communication device are implemented in the application software of the proc- 
essor block 17. Thus, a display 1 9 can also be used, in a way known as such, e.g. to display call information. The user 
can use a keypad 20 to control the positioning receiver and the wireless communication device. A codec 21 is used 
to code and decode audio signals. Furthermore, Fig. 6 shows a radio part 23 and a second antenna 23 for the wireless 
50 communication device. 

[0053] The present invention is not limited to the above-presented embodiments, but it can be modified within the 
scope of the appended claims. 

55 Claims 

1. A method for synchronizing a receiver (1) with a transmitted code-modulated spread spectrum signal, in which 
method at least one reference code (r(x)) is used : which corresponds to a code used in the modulation, and the 
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frequency shift of the transmitted signal and the code phase of the code used in the modulation are determined, 
characterized in that a correlation step is taken to form a correlation function matrix on the basis of the received 
signal and said reference code, a non-coherent search matrix is formed of said correlation function matrix, and 
that the method comprises at least one of the following steps: 

5 

modifying the elements of said non-coherent search matrix on the basis of at least one statistical property of 
the elements of said non -coherent search matrix, 

performing high-pass filtering before forming said correlation function matrix. 

w 2. The method according to claim 1 , characterized in that said non-coherent search matrix is two-dimensional, 
wherein one dimension is the frequency and the other dimension is the code phase of the reference code. 

3. The method according to claim 2, characterized in that in the modification of the elements, a parameter describing 
the deviation of signals is calculated from the values of the different code phases for each frequency value in the 

15 non-coherent search matrix, and said parameter describing the deviation of signals, multiplied with a scaling factor, 

is subtracted from each value of the code phase of said frequency value. 

4. The method according to claim 2 or 3, characterized In that in the modification of the elements, the average is 
calculated from the values of the different code phases of each frequency value in the non-coherent search matrix, 

20 and said calculated average is subtracted from each value of the code phase of said frequency value. 

5. The method according to any of the claims 1 to 4, characterized in that said correlation function matrix is formed, 
that after the formation of said correlation function matrix, the signal is sampled at a predetermined sampling 
frequency, and that the frequency of at least one spurious signal is determined in said signal after the sampling, 

25 to suppress the effect of the interference, which is at least one of the following: 

a spurious signal close to the frequency of the signal to be received, 
DC shift interference. 

30 6. The method according to claim 5, characterized in that the received signal is subjected to frequency conversion 
by mixing the received signal with a local oscillator signal, and that said determination of the frequency of at least 
one spurious signal is made on the basis of the frequency of said local oscillator signal. 

7. A system comprising a receiver (1) with means (6) for acquisition of a code-modulated spread spectrum signal, 
35 means (16) for generating at least one reference code (r(x)) corresponding to the code used in the modulation, 

and means (15) for determining the frequency shift of the transmitted signal and the code phase of the code used 
in the modulation, characterized in that the system also comprises correlation means (102-105) for forming a 
correlation function matrix on the basis of the received signal and said reference code, means (1 07-1 09) for forming 
a non-coherent search matrix from said correlation function matrix, and at least one of the following means: 

40 

means ( 1 5) for modifying the elements of said non-coherent search matrix on the basis of at least one statistical 
property of the elements of said non-coherent search matrix, 

means for performing high-pass filtering before the formation of said correlation function matrix. 

45 8. An electronic device comprising a receiver (1 ) with means (6) for acquisition of a code-modulated spread spectrum 
signal, means (16) for generating at least one reference code (r(x)) corresponding to the code used in the modu- 
lation, and means (15) for determining the frequency shift of the transmitted signal and the code phase of the code 
used in the modulation, characterized in that the electronic device (24) also comprises correlation means 
(102-105) for forming a correlation function matrix on the basis of the received signal and said reference code, 

so means (1 07-1 09) for forming a non-coherent search matrix from said correlation function matrix, and at least one 

of the following means: 

means ( 1 5) for modifying the elements of said non-coherent search matrix on the basis of at least one statistical 
property of the elements of said non-coherent search matrix, 
55 - means for performing high-pass filtering before the formation of said correlation function matrix. 

9. The electronic device (24) according to claim 8, characterized in that said non-coherent search matrix is two- 
dimensional, wherein one dimension is the frequency and the other dimension is the code phase of the reference 
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code. 

10. The electronic device (24) according to claim 9, characterized in that the means for modifying the elements in 
said non-coherent search matrix comprise means for computing the standard deviation from the values of the 

5 different code phases for each frequency value in the non-coherent search matrix, and means for subtracting said 

computed standard deviation, multiplied by a scaling factor, from each code phase value for said frequency value. 

1 1 . The electronic device (24) according to claim 9 or 1 0, characterized in that the means for modifying the elements 
in said non-coherent search matrix comprise means for computing the average from the values of the different 

10 code phases for each frequency value in the non-coherent search matrix, and means for subtracting said computed 

average from each code phase value for said frequency value. 

12. The electronic device (24) according to claim 9 or 10, characterized in that it comprises means for sampling the 
signal after the forming of said correlation function matrix. 

15 

13. The electronic device (24) according to any of the claims 8 to 12, characterized in that it comprises means for 
performing functions of a mobile station. 

14. The electronic device (24) according to any of the claims 8 to 1 3, characterized in that the receiver (1 ) is a satellite 
20 positioning receiver. 

15. The electronic device (24) according to any of the claims 8 to 14, characterized in that the correlation means 
(1 02-1 05) comprise a matched filter. 
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(54) A method, a system and an electronic device for synchronizing a receiver 



(57) The invention relates to a method for synchro- 
nizing a receiver (1) with a transmitted code-modulated 
spread spectrum signal. In the method, at least one ref- 
erence code (r(x)) is used, which corresponds to a code 
used in the modulation, and the frequency shift of the 
transmitted signal and the code phase of the code used 
in the modulation are determined. In the method, a cor- 
relation step is taken to form a correlation function matrix 
on the basis of the received signal and said reference 



code. The step is reiterated a plurality of times. The it- 
erations results are added non-coherently to form a non- 
coherent search matrix. Elements in said non-coherent 
search matrix are modified on the basis of at least one 
statistical property of the elements in said non-coherent 
search matrix, and/or high-pass filtering is performed 
before the formation of said correlation function matrix. 
The invention also relates to a system and an electronic 
device in which the method is applied. 
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- a correlation matrix is formed from the results of the 
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- a non-coherent search matrix is formed from the 
correlation matrix and 

- the elements of the non-coherent search matrix are 
modified according to at least one statistical property. 
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Method , system and electronic device for synchronizing a 
spread spectrum receiver through determination of frequency 
shift and spreading code phase, wherein 

- a correlation matrix is formed from the results of the 
correlation between the code and the received signal; 

- a non-coherent search matrix is formed from the 
correlation matrix and 

- highpass filtering is performed before forming the 
correlation matrix. 
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